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f"^ ' Abstract. The notion of a compact object immune to the horizon problem and 

comprising an anisotropic inhomogeneous fluid with a specific radial pressure behavior, 
i.e. the gravastar, is extended by introducing an electrically charged component. 
Einstein Maxwell field equations are solved in the asymptotically de Sitter interior 
O -1 ! where a source of the electric field is coupled to the fluid energy density. Two different 

| ' solutions which satisfy the dominant energy condition are given: one is the (5-shell 

model for which the analysis is carried out within Israel's thin shell formalism, the 
other approach — the continuous profile model — is solved numerically and the interior 
■ solutions have been (smoothly) joined with the Reissner-Nordstrom exterior. The 

effect of electric charge is considered, and the equation of state, the speed of sound 
and the surface rcdshift are calculated for both models. 
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1. Introduction 

^ ■ The search for solutions of Einstein's equations with different inputs for the physical 
content (dust, perfect fluid, etc.) and features (charge, anisotropy, viscosity, etc.) has 
been one of the steppingstones on the way to achieving comprehensible picture of today's 
Universe. Among solutions found so far one eminent position certainly belongs to black 
hole solutions with intriguing properties and characteristics (see for instance Ref. [1] 
for a recent review). However, some of the black hole features have been shown to 
be unattractive, apart from the problem that irrefutable observational evidence of the 
black hole existence has not yet been found. These reasons had inspired the search for 
alternative configurations which led to a solution dubbed gravastar, the gravitational 
vacuum star, by Mazur and Mottola [21 [3] . These spherically symmetric static solutions 
to the Einstein equations — candidates for highly compact astrophysical objects and in 
this sense alternatives to black holes — evolve from the segment of the de Sitter geometry 
in the center with the equation of state p + p = 0, proceed through a thin vacuum 
phase transition layer, avoid the event horizon formation, and swiftly match the exterior 
Schwarzschild spacetime. (Objects with asymptotically de Sitter core were considered 
also in earlier literature, see e.g. Refs [H El IB]-) One development of the gravastar 
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idea went in the direction of generalizing the interior and the exterior geometries and 
modeling the phase transition layer by a 5-shell [7J El El [101 EH E2l E3] . Another direction 
of development focused on replacing the 5-shell by a continuous transition of the energy 
density and pressure profiles from the asymptotically de Sitter center to the exterior 
Schwarzschild metric [HI [15] . The common feature of all realizations of gravastars 
is the anisotropy of principal pressures [161 EH EB]- Several astrophysically relevant 
aspects of gravastar solutions such as thermodynamic properties [10] . modes of quasi- 
normal oscillations [20] . and ergoregion instability [21] |22] were also discussed in the 
literature. 

Charged anisotropic models were considered in a number of papers [231 EU 1251 126] . 
Requiring T r r = for a charged anisotropic fluid, the solutions may represent 
electromagnetic mass models of electrically neutral systems |27J, since the total charge 
must be zero if one wants to have a configuration with the vanishing radial pressure at 
the surface. Another charged object with neutral core and the electric charge distributed 
on a 5-shell was investigated in Ref. [28]. In the same context the anisotropic charged 
spheres with conformal motion/symmetry were considered in Refs [251 [29]. Anisotropic 
charged spheres with (varying) cosmological constant were considered in Ref. [3D] . 

In this paper the gravastar picture is extended to include the effect of the electric 
charge as a natural step in gravastar investigations. Although astrophysical objects are 
essentially neutral, the problem of the electric charge in the phenomenological context 
could occur in the (strange) quark stars considerations, or accreting objects. Also, 
the influence of the electric charge on the spacetime curvature and other features of 
the Einstein-Maxwell system, could be seen in the context of the model of a classical 
charged massive particle. Therefore the electric charge extension in the gravastar context 
could be understood as a natural step in investigations. The solutions for the charged 
gravastar obtained here satisfy the dominant energy condition (DEC) everywhere and 
possess no horizons. The DEC requires the energy density to be non-negative, and the 
absolute values of the pressures not to exceed the value of the energy density. It is taken 
as the principal criterion for the viability of solutions. (The less stringent weak and null 
energy conditions are automatically satisfied if the DEC is satisfied [31].) The strong 
energy condition (SEC), which requires that the sum of the energy density and all of 
the pressures is non-negative (and guarantees the attractive character of gravity), is by 
definition violated in the de Sitter core of the gravastar [71 [12]. The interior metric is 
obtained by solving the Einstein-Maxwell equations in which the matter — anisotropic 
inhomogeneous charged fluid — through the Gauss law constraint of electrostatics serves 
as a source for the electric field. Such charged objects induce the Reissner-Nordstrom 
(RN) exterior metric which smoothly joins the regular interior solution. 

Two different charged gravastar models are presented. In Sec. [2] we consider a 
solution with a 5-shell charge distribution. We derive bounds on the configuration 
parameters that follow from the requirement that the DEC is satisfied and that the 
emergent equation of state leads to subluminal speed of sound. Physically viable 
solutions with unbounded surface redshift are shown to exist for near-to-extremally 
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charged <5-shell gravastars. 

Another solution presented in Sec. [3] stems from the notion of the anisotropic fluid 
with continuous pressure profiles and with an electrically charged component coupled 
to the fluid energy density [27J [32j [33J, [26j [M] . This choice simplifies the calculations, 
but does not restrict the overall generality since the 5-shell model of Sec. [2] can be 
understood as the model of perfectly conducting sphere with the electric charge expelled 
to the surface. The numerically obtained solutions reveal a strong dependence on the 
amount of the electric charge, but not on the particular shape of the assumed fluid 
energy density profile (as long as it obeys usual requirements [35]). For both solutions 
with different charge content the redshift is calculated, the equation of state (EOS) is 
inferred in the atmosphere of the gravastar where standard physics is expected to be 
valid [Hj. Also in this region the speed of sound is calculated to test the viability of 
results. Conclusions are given in Sec. HI 

2. The 5-shell model 

2.1. The model 

In the single 5-shell picture of the electrically charged gravastar, the interior geometry 
is a segment of the de Sitter (dS) spacetime, while the exterior is the asymptotically 
flat segment of the Reissner-Nordstrom (RN) spacetime. The two metrics are joined 
at a spherical timelike hypersurface representing the 'vacuum phase transition' layer 
of the original Mazur-Mottola model. Israel's thin shell formalism [36J formulates the 
conditions for the smooth joining of the metrics at a hypersurface. The first junction 
condition requires that the metric on the hypersurface induced by the metrics on its two 
sides be the same. The second junction condition requires that the extrinsic curvature 
of the hypersurface be the same when computed from the metrics on the two sides. 
If the second junction condition is not satisfied, then the joining is not smooth, but 
it is still allowed by the Einstein equations provided that a particular distribution of 
energy-momentum is introduced on the hypersurface. Such distribution is understood 
as 5-shaped in the complete spacetime. 

The dS metric of the gravastar interior and the RN metric of the exterior spacetime 
can both be written using the geometrized units and the coordinates x a = (t, r, ip) as 



where fi(r) is the so-called compactness function, k is the time coordinate scaling 
constant and dfl 2 is the metric on the unit sphere. The compactness vanishes at the 
regular center of a spherical body, and also at the asymptotically flat spatial infinity, 
while for all r it must remain less than unity in order to avoid formation of horizons in 
the spacetime. The junction hypersurface is specified with its radius which we denote 
r = a. The metric on the junction hypersurface induced from the general metric (0Q) can 
be written adopting the coordinates y a = (t, <p) as 




(1) 




(2) 
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The first junction condition is automatically satisfied provided that k dS (l — fids( a )) — 
&rn(1 ~ A i RN(fl)), where we choose to put fc RN = 1 as usual. The non- vanishing 
components of the extrinsic curvature of the junction hypersurface, computed from 
the general metric (pQ), are 



H'(a) 



* a 



//(a) 



(3) 



2^1 -//(a) 

(prime denotes derivative with respect to r) . According to the thin shell formalism, the 
hypersurface energy-momentum tensor is given by 

1 



St = - 



S7T 



(4) 



where K = K® and the square brackets indicate the discontinuity of the quantity across 
the hypersurface, in our case [/] = lim e ^ f(a + e) — f(a — e) = /rn(o) — /ds( a )- The 
non-vanishing components of the hypersurface energy-momentum tensor (j3J) computed 
from the extrinsic curvature ([3]) are interpreted as the surface energy density of the 
gravastar shell a and the isotropic surface tension 9, and are given by 



a= - Si 



1 - n(a) 



An a 



9 



QV 

0,a 



1 — /i(a) — a/x'(a)/2 



8n a Jl — /i(a) 



(5) 



(6) 



(note that by definition surface tension has the opposite sign of surface pressure). The 
compactness functions for the dS and the RN spacetimes can be written as 

Svrpo 2 2M Q 2 

Vds{r) = ^—r , //RN(r) = -, (7) 

6 r r z 

where p is the constant (dark) energy density of the dS spacetime and M and Q are 

the mass and the charge parameters of the RN spacetime. It is convenient to use the 

configuration variables defined as follows, 

K = \Q/M\ , x ee /i dS (a) G (0, 1) , 1/ ee /i RN (a) G (0, 1) . (8) 

(The definition of x and y as the compactness at the interior and at the exterior side of 
the gravastar shell is compatible with the notation used in Ref. |12j.) For the charge to 
mass ratio k < 1, the RN spacetime has the horizons at r± 

as k 



M{1 ± VI - k 2 ) which, 

1, merge into a single (extremal) horizon of the ERN spacetime. For k < 1, 
in order to exclude the horizons from the geometry, the gravastar surface radius must 
satisfy a > r A 



M(l + VT 



K 2 ) 



For k > 1, the compactness is bounded from above 



with its maximum k~ 2 < 1 at r = Mk 2 and there are no horizons in the RN spacetime. 
However, this upper bound also forbids arbitrarily compact configurations so the k > 1 
configurations will not be further considered. Assuming k < 1 for the rest of this section, 
and using the configuration variables x, y, the surface energy density o and the surface 
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tension 9 of the gravastar shell are given by 

Q = i-2x + i - k 2 - yi - ^ 

87ra \/l — x 8ir a k 2 y/1 — y 
2.2. Energy conditions 

We now want to determine the region in the x, y configuration space where the DEC 
(and therefore also the NEC and the WEC) is satisfied for different values of k (as 
mentioned earlier we do not require the SEC to be satisfied). The DEC requires the 
surface energy density to be non-negative and the absolute value of the surface tension 
(or pressure) not to exceed the value of the surface energy density, i.e. 

< \6\ < a . (11) 

The condition a > with the surface energy density ([9]) implies x < y, i.e. the 
compactness of the interior metric must be less than the compactness of the exterior 
metric at the shell. Therefore the configurations in lower right triangle of the x, y plane, 
shown in Fig. [1] (left plot), are not allowed by the DEC. Inspection of ffTU]) reveals that 
9 < for all values of x, y G (0, 1) and k G (0, 1], meaning that in all configurations we 
are considering the gravastar shell is under (positive) surface pressure. The condition 
\9\ < o imposed by the DEC therefore reduces to 9 + o > 0, which for (TTUT) further 
reduces to an inequality of the form 

x<f(y-,K) , ye (0,1), (12) 

where the function f(y; k) can be obtained analytically, but is a rather complicated 
expression so we discuss only its most important features. In the case k = 0, 
corresponding to the electrically neutral gravastar, (fl"2l can be written as 



60 - 3Qy - 25y 2 - (6 - SyWlOO - 2/(124 - 25y) 
^ /(2/;0) = 128(1^ • (13) 

The function f(y; 0) has nodes at y — and y = 24/25, for y e (0,24/25) it satisfies 
< f{y; 0) < y with a single maximum at y = 4/5, /(4/5; 0) = (19- v / 105)/32) ~ 0.274, 
see the innermost shaded 'D-shape' in Fig. [1] (configurations within the 'D-shape' 
are allowed by the DEC). This means that the highest surface compactness on the 
exterior side of the shell of the electrically neutral gravastar satisfying the DEC is 
y = A*rn(o) — > 24/25 with x = /ids (a) ~" > 0, i.e. with vanishing (dark) energy density 
in the interior. This also means that, if the DEC is to be satisfied, increasing the 
(dark) energy density in the interior does not lead to higher surface compactness of the 
electrically neutral gravastar. 

For the charge-to-mass ratio k G (0,1), the DEC satisfying 'D-shaped' region in 
the x, y parameter space is larger than that corresponding to the neutral configuration, 
but the boundary f(y; k) is qualitatively similar to the k = case. As shown in the 
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Fig. [T] (left plot), the upper y-node of f(y; k) (the 'D-shape') for k G (0, 1) remains less 
than unity, implying that the surface compactness of the gravastar satisfying the DEC 
is bounded from above below unity. Again, an increase of the (dark) energy density 
in the gravastar interior corresponds to a decrease of the upper bound on the surface 
compactness. 

A qualitatively different situation is reached at K — 1. Here we obtain 
x < f(y; 1) = ^ ( 19 - 4^1^ + 4y 

- ^121 + 104 V /T^ - 8y (l5 + - 2yj j , (14) 

see Fig. [T] (left plot). The boundary of the DEC-allowed region is no longer a 'D-shape'. 
The upper bound on the surface compactness is no longer present, which means that 
the surface compactness can reach unity arbitrarily close without violating the DEC. 
The gravitational redshift of light emitted from the surface of such a body could be 
arbitrarily high. 



2.3. Equation of state and causality 

In a perfect fluid with the equation of state p = p(p), where p is the pressure and p 
is the energy density, the quantity c 2 = dp/ dp can be interpreted as the speed of the 
propagation of sound (squared). Negative c 2 is usually interpreted as an indication 
of instability, or at least as the impossibility of propagation of sound, while c 2 > 1 
indicates superluminal speed of sound or violation of causality. In the context of the 5- 
shell gravastar model, the matter comprising the shell can be understood as the perfect 
fluid in 2+1 dimensions. In Ref. [7] a procedure has been developed that allows one to 
extract the equation of state 9 = 9(a). Starting from the solutions obtained above one 
can extract the equation of state 9 = 9(a) from which the speed of sound (squared), 

2 d9 d9/dr 



da da/dr 



(15) 



becomes 



(k 4 - 3k 2 + 2 + ((3 - V )k 2 - 2) VI - tfy) 



c 2 = \=i x - y \ _: - (i6) 

2/? 4 VT^ - 2k 2 v / T^ (k 2 (1 - 2y) + 1 - VI - ^y) 



As shown in Fig. [T] (right plot), for k G [0, 1) we have the c 2 < (non-propagating) 
regime in the lower-right region of the x, y square, the c 2 > 1 (superluminal) regime is in 
the upper- left, while the physically acceptable region is confined in a 'bended triangular 
region' bounded by the solid curve from above (c 2 = 1), and the corresponding dashed 
curve from below (c 2 = 0). Starting from the lowest pair of curves (dashed and solid) 
bounding the shaded c 2 G [0, 1] (allowed) region for the electrically neutral gravastar, as 
k increases we see that the 'bended triangle' shifts toward higher surface compactness y, 
but also becomes narrower. In the limit «;—>!, the c 2 G [0, 1] (allowed) region shrinks 




Figure 1. The (5-shell charged gravastar configurations shown in the x, y plane. Energy 
conditions (left plot): The configurations with x > y (lower-right triangle below the 
dashed line) would require negative surface energy density a and are thus ruled out 
by the DEC. In the 'D-shaped' regions bounded by solid lines the condition |0| < a 
required by the the DEC is satisfied. Going from left to right, the solid lines correspond 
to k = \Q/M\ = (shaded region), and subsequent lines to n = |, |, || , |p The 'open 
top' solid line corresponds to n — 1 (ERN). Bullets indicate highest x that can be 
obtained with certain k. Causality of EOS (right plot): Pairs consisting of a dashed 
and a solid line bound the regions in which cl? € [0, 1] for different values of k — \Q/M\. 
In the regions below the dashed lines Cg < 0, and above the solid lines Cg > 1. The 
lowest pair of dashed-solid lines bounding the shaded 'bended triangle' corresponds 
to k = 0, and subsequent pairs of lines correspond to to k = |, |, || , |j. For n = 1 
(ERN) only the dashed line is shown, since for all configurations ci? < 0. 



so that the < (non-propagating) region takes over the whole x, y square. Therefore, 
the k — 1 configurations (ERN), do not satisfy the Cg G (0, 1) requirement. 

3. Continuous profile models 

3.1. The model 

Here we take the similar approach but in the Einstein-Maxwell system governed by the 
equations 

Rl - = SttT?, = A-Kf , (17) 

where Rp is the Ricci tensor, R = R®, = T( fluid )^ + T( char g e )^ is the energy-momentum 
tensor of the anisotropic fluid and the electrostatic field. The electromagnetic part of 
the energy- momentum tensor is given by T( charge )^ = -^(F^F^ — ^5pF al3 F a/3 ), F a/3 is 
the electromagnetic field strength tensor and j a the source four-current. We use the 
coordinates x a = (t, r, t?, <p) as before, and write now the general spherically symmetric 
static metric simply as 

ds 2 = g tt (r) dt 2 + g rr .(r) dr 2 + r 2 dtt 2 . (18) 
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Assuming A a = S^A^r) for the electromagnetic four-potential, the only non- vanishing 
component of the field strength tensor F af3 = d a A@ — dpA a is F tr = —F rt = —d r A t , and 
the locally measured radial electric field is E = d r A t /y/—gu g rr - The four-velocity 
of the static fluid element is u a = 5f /y/— gtt, and the four-current of the electric 
charge with density cr(r) is j a = a u a = 5fo~/^—gu- The Maxwell equation reads 
(r 2 E(r))' = An r 2 cr/y/g^, (prime denoting the derivative with respect to r). Integrating 
from the center to some radius r one obtains 

E(r) = — 5- , where q(r) = / Anr' 2 cr(r') \/g rr (r')dr' (19) 
r 2 Jo v 

is the amount of the electric charge within the sphere of radius r. The electromagnetic 
part of the energy momentum tensor is 

T( charge )^(r) = ^ diag(- -, +, +) , (20) 

and the part due to the anisotropic fluid is 

T( flu id)^(r) = diag (-p(r), p\\(r), p ± (r), pj_(r)) , (21) 

where p is the energy density, p\\ is the radial and p± is the transverse pressure. Writing 
the metric components gu and g rr in terms of two functions, m and as 

g tt (r) = - (1 - 2m(r)/r) e 2 * (r) , g rr (r) = (1 - 2 m (r)/r)^ , (22) 

the Einstein equations give three equations: 

m'(r)=A7rr 2 p(r) + ^£- , (23) 

2r z 

(p(r) +p||(r))(m(r) + 47rr 3 p||(r) — g 2 (r)/2r) 
P||( r ) = r 2 (l -2m(r)/r) 

+ - i {r) q{T ) = ; + - (P±(r) - P|,(r)) , (24) 
rvl — 2m(r)/r r 

nr)= ^ r M +p\\<r) , (25) 

1 — 2m(r)/r 

Eq. (|23|) defines the 'mass function' m(r) and (1241) is the version of the well-known TOV 
equation for the anisotropic and electrically charged fluid. 

The electrically neutral gravastar solutions with continuous pressure profiles 
satisfying the general requirements of Ref. [13] were constructed in Ref. [15] from the 
above system (with q = 0) together with these essential ingredients: 

• continuous fluid energy density profile p(r) outwardly non-increasing and vanishing 
at the surface radius r = R, for example 

p(r) = po (1 - (r/R) n ) ; n>2, (26) 

where po is the central fluid energy density, 

• fluid pressure anisotropy Ansatz of the form 

!! ^-U*P) t &!}'■• kJ> - u (27) 

where a is the measure of anisotropy strength, 



Electrically charged gravastar configurations 



9 



• behavior of the pressures - boundary conditions, 

P||(0)=p_L(0) = -p(0) , m (R)=p ± {R) = 0. (28) 

For the electrically charged model, we adopt the above ingredients, and couple two 
density profiles in the following way 

a(r) = ep{r) - 2m(r)/r , (29) 

with e constant [33 [23J [23 133] • The equations (123 H25j) together with the Maxwell 
equation that can be written as 

q'(r) = 4vr r 2 a{r) 1 = , (30) 

Jl — 2m(r)/r 

close the system in four unknowns: the two metric components m(r) and \I/(r) , the fluid 
radial pressure p\\(r) and the charge q(r). 

At r = i? we match the interior solution to the exterior segment of the Reissner- 
Nordstrom spacetime with the metric 

^ = -( 1 -™ + $l)dt 2 +(l- — + %Y 1 + *t , (31) 
V r r 2 J \ r r 2 J 

where M and Q are the mass and the charge parameters of the RN spacetime. 

The interior mass function can always be written as the sum of two functions, 

m(r) = m p (r) + m q (r) , (32) 

where 

(\ f r A ,2 ( s r q 2 {r') , foo\ 

m p\ r ) = = / 47rr p[r)dr , mJr) = / — ar , (33) 



2r' 2 

subscripts suggesting the origin of contributions: fluid and charge. The next step in the 
solution process is to solve the TOV ff24l) for the radial pressure. (We note here that 
the equations are invariant with respect to the change of sign of the electric charge.) 
Acceptable solutions of the equations, for a given combination of k, I, n in (126]) and (J27j) 
depend on a po and e. A step in cornering their values is connected with the smooth 
joining of the metrics at the gravastar surface 

O 2 

m(R) = nrip(R) + m q (R) = M — and q(R) = em p (R) = Q . (34) 

2R 

Further denoting k = \Q/M\ as in Sec. [2} the above conditions combine to give a 
quadratic in e 

from which we get the upper bound on the central fluid energy density p . The bound 
is saturated at e — £x,2 = 2k. Two values of e lead to the desired charge-to-mass ratio 
k = \Q/M\, one of which leads to solutions violating the DEC. 
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Figure 2. Continuous profile gravastar solution inducing the ERN (re = \Q/M\ = 1) 
spacetime: energy density ([25)1 with n = 2 and pressure anisotropy (f2"T|) with A; = / = 1 . 
Using s = ei = e-2 = 2k gives poR 2 = 21/327T, M = Q = 7R/10, surface compactness 
fi = 91/100, surface redshift = 7/3, and pressure anisotropy strength a 2 /12 ~ 0.69. 

3.2. Results 

The solution displayed in Fig. [2] corresponds to the maximum allowed (fluid) density 
Pq = 21/ \32rrR 2 ). The presented solution corresponds to the ERN case and it is 
evident that the radial and transverse pressures obey the DEC while the compactness is 
safely protected from reaching unity. The gravastar configurations obeying the energy 
conditions [12] show no quasi black hole behavior as described in Ref. [38] . 

In the context of compact objects an important source of information and 
classification criterion is the redshift, which in the uncharged gravastar version was 
shown to be within expectations for compact objects [15J. The surface redshift given by 

Z^ = -l + \g tt {R)\- 1 l\ (36) 

for the perfect fluid spheres is less or equal to 2. The anisotropy alters this value and 
provides a significant increase up to 3.84 [391 HOI HI], while the uncharged gravastar 
solutions with the continuous pressure [15] give values between 1.23 and 1.71. The 5- 
shell model for the ERN case allows the unbounded surface redshift value, whereas the 
continuous density model gives the surface redshifts that are larger than their uncharged 
counterparts, e.g. in the ERN case shown in Fig. [2] the surface redshift is ~ 2.33 
with maximum compactness ~ 0.92. 

The results quoted so far do not depend qualitatively on the particular choice of the 
profile function for the fluid/charge density distribution: for n — 4, the corresponding 
maximal central density is po — 585/(12167ri? 2 ), and the maximal (internal) compactness 
here is /i = 0.93, so we get basically the same form of solutions. 

For the general RN solutions the r + horizon is safely positioned within the gravastar 
surface: the interior compactness fi(r) = 2m(r)/r = 2(m p (r) +m q (r)) /r has to be always 
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Figure 3. The behavior of fluid radial pressure in the gravastar atmosphere and crust 
for energy density (|26|) with n — 2 and pressure anisotropy (|27|) with k = I = 1, 
for a sequence of solutions with poi? 2 = 21/327T and n = \Q/M\ = 0,0.5,0.75,0.99: 
Pressures are shown as functions ofr/R (upper plot) and as functions of the fluid energy 
density (lower plot, solid lines). In the polytropic equation of state approximating the 
behavior of the pressures in the outermost layer of the atmosphere, p = kp 1+1 / np , the 
polytropic index is n p ~ 1 for all solutions and k decreases with n (thin dashed lines). 



smaller then unity. This gives 2m(R) < R and since the RN mass parameter M is given 
by (1321) then M — > m(R) for Q — > giving r + — > 2M for Q — > 0. Therefore 

r + -> 2M -> 2m(i2) < i2 or r + < R , (37) 

and we have no "naked horizon" for the charged gravastar. 

As a result of our approach to the charged gravastar problem we are in the position 
to correlate the pressure and the fluid density. The correlation changes from the usual 
matter at low densities to that associated with de Sitter geometry at high densities near 
the center of configuration, and this change is well behaved. It is important to note 
that in the atmosphere of the gravastar (see [HI US]) the radial pressure has a negative 
gradient, so in this region a star-like behavior is expected. This behavior is shown in 
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Fig. [21 In the upper panel it is interesting to note the change of the behavior of the 
radial pressure due to the presence of charge: by increasing the charge the fluid pressure 
decreases. Since the electric energy density adds on to the fluid density, the total mass 
and the compactness are increasing as well. 

In the lower panel of Fig. [3] several solutions are presented all of which show the 
polytropic pattern in the low density region, 

p(p) = kp 1 = kp 1+1 ^. (38) 

For all solutions the polytropic index is n p ~ 1 and k decreases with the charge-to-mass 
ratio k. 

The speed of sound calculated from the EOS is always subluminal for all solutions 
in Figs [2] and [3J and is less than in the case of the uncharged gravastar. Recall that the 
negative speed of sound squared excludes the ERN case in the S— shell model. 

4. Conclusions 

In this paper two versions of the charged gravastar solutions have been presented. 
Within the 5-shell model, solutions that satisfy the DEC have been obtained and 
it has been shown that unbounded surface redshift is possible only in the solutions 
inducing the ERN spacetime. However, the analysis of the speed of sound computed 
from the resulting equation of state revealed that the ERN solutions are excluded. 
In the continuous profile approach, solutions are determined by values of the central 
fluid energy density p and the charge-to-mass ratio k. Close to the surface where the 
pressures are non-negative the equation of state is constructed. It shows the polytropic 
behavior with the index n p ~ 1. The speed of sound is subluminal. 

The obtained results help further understanding of the gravastar concept, either as 
the possible alternative to the black hole and phenomenological cosmological objects, 
or simply as interesting solutions found within the theory of General Relativity. 

Acknowledgments 

We thank the referees for suggestions that helped us to improve the paper. We thank 
Andrew DeBenedictis for reading the manuscript and for useful remarks. This work has 
been supported by the Croatian Ministry of Science, Education and Sport under the 
project No. 036-0982930-3144. 

References 

[1] A. DeBenedictis, "Developments in Black Hole Research: Classical, Semi-Classical, and 

Quantum", in M. N. Christiansen and T. K. Rasmussen, eds., Classical and Quantum Gravity 
Research. Nova Sci. Pub., 2007. |0711.2279 [gr-qc]| 

[2] P. O. Mazur and E. Mottola, "Gravitational condensate stars," gr-qc/0109035. 

[3] P. O. Mazur and E. Mottola, "Gravitational vacuum condensate stars," Proc. Nat. Acad. Sci. 
Ill (2004) 9545, ! gr-qc/0407075[ 



Electrically charged gravastar configurations 



13 



E. B. Gliner, "Algebraic properties of the energy- momentum tensor and vacuum-like states of 
matter," Sov. Phys. JETP 22 (1966) 378. 

A. D. Sakharov, "The Initial Stage of an Expanding Universe and the Appearance of a 
Nonuniform Distribution of Matter," Sov. Phys. JETP 22 (1966) 241. 

I. Dymnikova, "Vacuum Nonsingular Black Hole," Gen. Rel. Grav. 24 (1992) 235. 
M. Visser and D. L. Wiltshire, "Stable gravastars - an alternative to black holes?," Class. 
Quantum Grav. 21 (2004) 1135, |gr-qc/0310107| 

B. M. N. Carter, "Stable gravastars with generalised exteriors," Class. Quantum Grav. 22 (2005) 
4551, gr-qc/0509087[ 

N. Bilic, G. B. Tupper, and R. D. Viollier, "Born-Infeld phantom gravastars," J. Cosmol. 
Astropart. Phys. 02 (2006) 013, |astro-ph/0503427[ 

F. S. N. Lobo, "Stable dark energy stars," Class. Quantum Grav. 23 (2006) 1525, 
|gr-qc/0508115| 

F. S. N. Lobo and A. V. B. Arellano, "Gravastars supported by nonlinear electrodynamics," 



Class. Quantum Grav. 24 (2007) 1069, gr-qc/0611083 



D. Horvat and S. Ilijic, "Gravastar energy conditions revisited," Class. Quantum Grav. 24 (2007) 

5637, 0707.1636 [gr-qc] [ 
P. Rocha, A. Y. Miguelote, R. Chan, M. F. da Silva, N. O. Santos, and A. Wang, "Bounded 

excursion stable gravastars and black holes," J. Cosmol. Astropart. Phys. 06 (2008) 025, 

|0803.4200 [gr-qc]| 

C. Cattoen, T. Faber, and M. Visser, "Gravastars must have anisotropic pressures," Class. 



Quantum Grav. 22 (2005) 4189, gr-qc/0505137 



A. DeBenedictis, D. Horvat, S. Ilijic, S. Kloster, and K. S. Viswanathan, "Gravastar solutions 
with continuous pressures and equation of state," Class. Quantum Grav. 23 (2006) 2303, 
|gr-qc/0511097[ 

R. L. Bowers and E. P. T. Liang, "Anisotropic spheres in general relativity," Astrophys. J. 188 
(1974) 657. 

L. Herrera and N. O. Santos, "Local anisotropy in self-gravitating systems," Phys. Rep. 286 
(1997) 53. 

L. Herrera, J. Ospino, and A. di Prisco, "All static spherically symmetric anisotropic solutions of 

Einstein's equations," Phys. Rev. D 77 (2008) 027502, [0712.0713 [gr-qc][ 
A. E. Broderick and R. Narayan, "Where are all the gravastars? Limits upon the gravastar 



model from accreting black holes," Class. Quantum Grav. 24 (2007) 659, gr-qc/0701154 
C. B. M. H. Chirenti and L. Rezzolla, "How to tell a gravastar from a black hole," Class. 

Quantum Grav. 24 (2007) 4191, |0706 . 1513 [gr-qc]] 
V. Cardoso, P. Pani, M. Cadoni, and M. Cavaglia, "Ergoregion instability of ultracompact 

astrophysical objects," Phys. Rev. D 77 (2008) 124044, [0709.0532 [gr-qc] | 
C. B. M. H. Chirenti and L. Rezzolla, "On the ergoregion instability in rotating gravastars," 

Phys. Rev. D 78 (2008) 084011, [0808 . 4080 [gr^^c]] 
T. Singh, G. P. Singh, and A. M. Helmi, "New solutions for charged anisotropic fluid spheres in 

general relativity," N. Cim. B 110 (1995) 387. 
R. Roy, A. A. Rangwala, and N. C. Rana, "Charged anisotropic fluid spheres in general 

relativity," Indian J. Pure Appl. Math. 27 (1996) 1119. 
J. Krishna Rao, M. Annapurna, and M. M. Trivedi, "Static charged spheres with anisotropic 

pressure in general relativity," Pramana 54 (2000) 215. 
L. Herrera and V. Varela, "Finite self-energy of pointlike sources for the Reissner-Nordstrom 

metric," Gen. Rel. Grav. 28 (1996) 663. 
0. Gr0n, "A charged generalization of Florides' interior Schwarzschild solution," Gen. Rel. Grav. 

18 (1986) 591. 

V. Varela, "Neutral perfect fluids and charged thin shells with electromagnetic mass in general 
relativity," Gen. Rel. Grav. 39 (2007) 267, |gr-qc/0604108[ 



Electrically charged gravastar configurations 



14 



L. Herrera and J. Ponce de Leon, "Isotropic and anisotropic charged spheres admitting a 

one-parameter group of conformal motions," J. Math. Phys. 26 (1985) 2302. 
S. Ray and S. Bhadra, "Energy density in general relativity: A possible role of cosmological 



constant," Phys. Lett. A 322 (2004) 150, gr-qc/0212119 



M. Visser, Lorentzian Wormholes: From Einstein to Hawking. AIP Press, 1996. 

J. Ponce de Leon, "General relativistic electromagnetic mass models of neutral spherically 

symmetric systems.," Gen. Rel. Grav. 19 (1987) 797. 
J. Ponce de Leon, "Limiting configurations allowed by the energy conditions," Gen. Rel. Grav. 

25 (1993) 1123. 

J. Ponce de Leon, "Electromagnetic mass-models in general relativity reexamined," Gen. Rel. 

Grav. 36 (2004) 1451, |gr-qc/0310117| 
T. W. Baumgarte and A. D. Rendall, "Regularity of spherically symmetric static solutions of the 

Einstein equations," Class. Quantum Grav. 10 (1993) 327. 
W. Israel, "Singular hypersurfaces and thin shells in general relativity," Nuovo Cimento B 44 

(1966) 1. 

F. I. Cooperstock and V. de La Cruz, "Sources for the Reissner-Nordstrom metric," Gen. Rel. 
Grav. 9 (1978) 835. 

J. P. S. Lemos and O. B. Zaslavskii, "Black hole mimickers: regular versus singular behavior," 

Phys. Rev. D 78 (2008) 024040, [0806 . 0845 [gr-ljc]] 
B. V. Ivanov, "Maximum bounds on the surface redshift of anisotropic stars," Phys. Rev. D 65 

(2002) 104011, 1^7^0/02010901 
D. E. Barraco, V. H. Hamity, and R. J. Glciscr, "Anisotropic spheres in general relativity 

reexamined," Phys. Rev. D 67 (2003) 064003. 
S. Karmakar, S. Mukherjee, R. Sharma, and S. D. Maharaj, "The role of pressure anisotropy on 

the maximum mass of cold compact stars," Pramana 68 (2007) 881,|0708~3305 [gr-qc] [ 



